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ABSTRACT 



Aims. We present the results of flux density monitoring of PSR B0329+54 at the frequency of 4.8 GHz using the 32-meter TCfA 
radiotelescope. The observations were conducted between 2002 and 2005. The main goal of the project was to find interstellar 
scintillation (ISS) parameters for the pulsar at the frequency at which it was never studied in detail. To achieve this the 20 observing 
sessions consisted of 3-minute integrations which on average lasted 24 hours. This gave us sufficient sensitivity to all types of flux 
density variations over a wide range of timescales. The character of the observations makes our project unique amongst other ISS 
oriented observing programs, at least at high frequency. 

Methods. Flux density time series obtained for each session were analysed using structure functions. For some of the individual 
sessions as well as for the general average structure function we were able to identify two distinctive timescales present, the timescales 
of diffractive and refractive scintillations. To the best of our knowledge, this is the first case when both scintillation timescales, 
'diss = 42.7 minutes and tj^ss = 305 minutes, were observed simultaneously in a uniform data set and estimated using the same 
method. 

Results. The obtained values of the ISS parameters combined with the data found in the literature allowed us to study the frequency 
dependence of these parameters over a wide range of observing frequencies, which is crucial for understanding the ISM turbulence. 
We found that the Kolmogorov spectrum is not best suited for describing the density fluctuations of the ISM, and a power-law spectrum 
with/? = 4 seems to fit better with our results. We were also able to estimate the transition frequency (transition from strong to weak 
scintillation regimes) as Vj. = 10.1 GHz, much higher than was previously predicted. We were also able to estimate the strength of 
scattering parameter u = 2.67 and the Fresnel scale as rp = 6.7 x 10* meters. 

Key words, stars: pulsars: individual: PSR B0329-I-54 - ISM: structure - plasma - turbulence 
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1. Introduction 

Scintillations of the pulsar signal is a phenomenon that was 
identified shortly after the discovery of pulsars themselves 
(Hewish et al. I1968I I. From an observer point of view pulsars 
are point-like sources, and presumably are intrinsically stable 
sources of radiation. Their radiation, however, undergoes scatter- 
ing process in the interstellar plasma on its way to the observer 

Depending on the strength of the scattering that the radiation 
of a given pulsar undergoes, the interstellar scintillations (ISS) 
will result in different phenomena observed for the pulsar (see 
Lorimer & Kramer, '2005' and the references therein for a recent 
review). For a nearby pulsar, or when the radiation encounters 
only small amounts of the scattering material (and if the observ- 
ing frequency is higher than ca. 1 GHz), the disturbances of the 
wavefront will be minimal, resulting in weak scintillations. 

In case of strong scattering (a distant pulsar and/or low ob- 
serving frequency), the phase perturbations become large, which 
leads to a strong modulation of the pulsar signal. In this strong 
scintillation regime one has to take both diffraction and refrac- 
tion effects into account. The diffractive interstellar scintilla- 
tions arise from small (10^ H- 10** meter) inhomogeneities in the 
interstellar medium, which introduce observed flux density vari- 
ations with the timescale that can be estimated as fpiss ~ /'^ d^^'^ 



(assuming a Kolmogorov spectrum of the ISM turbulence, see 
Rickettil990j Romani et al. [T986i ). 
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The refractive interstellar scintillations (RISS; Sieber[T982| 
Rickett at al. 119841 ) arise from larger scale density fluctuations 
(10'" H- 10'-^ meters). As theories predict, the RISS timescale de- 
creases with observing frequency, roughly as fRiss ~ /^'^ d^'^- 

For any given pulsar the values of the mentioned parame- 
ters will hence change with the observational frequency. Starting 
from very low frequencies, we will have short but increasing 
DISS timescales, and very long but decreasing RISS timescales, 
which will be accompanied by increasing RISS modulation 
(DISS modulation will be always close to unity). The rate at 
which these parameters change can be used to decide between 
various theories describing the ISM turbulence spectrum, be- 
cause their frequency dependence power-law predictions are tied 
to the model used (see Romani et al l 19861 Bhat et al. |1999D . 

As the observing frequency increases, both scintillation 
timescales should be closer, and the refractive modulation in- 
dex should increase towards unity (which is best illustrated by 
Lorimer & Kramer. .2005. on their Figure 4.5b). At a certain fre- 
quency, called the transition frequency both timescales become 
equal, refractive scintillations reach maximum modulation and 
the character of the scintillations in general changes. This is the 
frequency at which the transition from the strong to weak scin- 
tillation regimes happens. 
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Most of the pulsar scintillation observations conducted so far 
were performed at the frequencies of 1 .5 GHz and below, there 
is only a handful of papers describing ISS observations at higher 
frequencies (Malofeev et al. 119961 Kramer et al. 119971 Shishov 
et al. I2003I I. The attempts to estimate the spectrum of the ISM 
turbulence via the measurement of ISS parameters were limited 
to the frequency range from ca. ~ 100 MHz to 1.5 GHz. This 
is both because pulsars are steep spectrum sources and are con- 
siderably weaker at higher frequencies, and the increasing DISS 
timescale which means one needs longer data sets to be able to 
perform measurements in a way similar to what is done at low 
frequencies. Both of these reasons result in a much longer ob- 
serving time to be able to obtain reliable measurements. Because 
there are only few radiotelesopes capable of high-frequency ob- 
servations, projects like this are very hard to implement. 

The way to at least partially overcome these obstacles is to 
limit the observations to very strong pulsars, which in turn low- 
ers the telescope-size requirements. One of the best candidates 
for these observations is definitely PSR B0329+54, the strongest 
radio pulsar in the northern hemisphere. Owing to its bright- 
ness, it was observed by many various ISS projects in the past 
(see the Discussion section in this paper, and Table [3] for a full 
list of references). Amongst them was a project conducted by 
Malofeev et al. ( 119961 ), which involved fairly short (ca. 130 min- 
utes) observations at 4.75 and 10.55 GHz, based on which the 
authors concluded that for PSR B0329+54 (with DM = 26.833) 
the transition frequency is ca. 3 GHz. This would mean that at 
higher frequencies this pulsar should be in the weak scintilla- 
tion regime, and this statement was contradicting our personal 
experience with that source. The flux density of PSR B0329+54 
measured for the purposes of estimating pulsar spectrum (see 
Maron et al. 2000 for an overview of that project) was showing 
very strong variations, and in individual measurements it ranged 
from a few to over a hundred milli-Janskys. 

To solve this problem we decided to perform a long-time 
monitoring project of PSR B0329+54 at the frequency of 
4.8 GHz. We were kindly given a chance of extensive use of 
Toruri Centre for Astronomy 32-meter radiotelescope (TCfA, 
located near Toruh, Poland), which is equipped with a cooled 
5 GHz receiver and the Penn State Pulsar Machine II. We used it 
to perform numerous long-duration observing sessions between 
2002 and 2005. This paper summarizes our results. 

2. Observations and analysis 

We have carried out a flux density monitoring of PSR B0329-I-54 
during three years using the 32-meter Toruh Centre for 
Astronomy radiotelescope. Observations were performed be- 
tween July 2002 and June 2005 at the frequency of 4.8 GHz. 
Twenty epochs of flux density monitoring at 4.8 GHz have been 
collected (see Table 1). Typical observations lasted about 12 - 30 
hours in 3 minute scans. During the first two years of our cam- 
paign each epoch was at least 24 hours; for the third year we 
observed this pulsar only occasionally, for about 12 hours each 
epoch. In addition, the first observing session lasted a full 5 days, 
and some other were almost 2 day long. 

The total intensity signal was obtained using the Penn State 
Pulsar Machine II (PSPM II). This pulsar back-end is basi- 
cally a 64 channel fast radio-spectrometer, which was designed 
to conduct timing observations and search for pulsars (see 
Konacki et al.'1999). The PSPM II is currently equipped with ad- 
ditional system that allows mean flux calibration. A dual-circular 
polarization signal from the total bandwidth of 192 MHz is fed 
into the backend, and then split into 64 channels (of 3 MHz width 



each, separately for each polarization). Signals from both po- 
larizations in each channel are then added, and sampled with a 
rate of 4096 samples per period, thanks to a system that allows 
dynamical changes of the sampling rate during the observation. 
This makes data much easier to integrate over time, which in fact 
is done in real-time. At the same time the noise diode signal is 
injected once per pulsar period, always on the same phase, which 
was at least 0.3 different from the actual pulse phase. As a result 
we have got 64 integrated pulse profiles with a calibration mark 
for each single observation. Data from all channels were dedis- 
persed and integrated over the whole bandwidth in an off-line 
procedure. 

To obtain flux densities, we carried out regular calibration 
measurements using the same signal of a noise diode, which was 
compared to the flux density of known continuum calibrators. 
Pulsar observations were interrupted regularly for the calibration 
procedure roughly every 6-8 hours (such calibration breaks took 
about one to one and a half hours), and sometimes there were 
unforeseen breaks owing to the telescope failures and hardware 
problems. 

2.1. Flux density measurements 

Figure [1] shows four of the total 20 flux time-series as observed 
during our campaign. Subplot (a) represents the flux variation 
during the longest observing session, i.e. the first epoch in our 
campaign, that lasted for five days. Evidently, flux variations 
happen at more than one timescale; the detailed analysis is pre- 
sented in the next sections. Subplot (b) shows the epoch with 
the highest measured flux value in our campaign. The maximum 
averaged flux for a single 3-minute integration reaches almost 
200 mJy, which means more than 10-times the average value 
{{F} = 17.4 mJy, see TablelT]i. This value is also close to the av- 
erage flux of PSR B0329-I-54 at 1400 MHz (203 mJy according 
to ATNF pulsar catalogue). An example of one of the worst qual- 
ity sessions is shown in subplot (c). There are many gaps in the 
data taken at that epoch, some are the regular calibration breaks, 
but others are caused bye the telescope/hardware failures. With 
all those gaps and some uncharacteristic flux variation it was 
not possible to conduct a proper analysis for this epoch (con- 
front Table[T]and Section l24l for detailed explanation). The last 
subplot (d) shows the flux variations for the penultimate epoch. 
Again, the flux variations are not typical for our data, the value 
stays well above the average for a very long time (more than 300 
minutes). This kind of variability makes some kinds of statistical 
analyses unreliable, i.e. we were not able to find the scintillation 
timescales because the structure function of the flux time series 
for this data set does not saturate. At the same time one has to 
note that this behaviour is quite typical for PSR B0329-H54, at 
least in our data. If one confronts the flux variations from sub- 
plot (d) with the longest time series (subplot (a)), one can easily 
find a matching behaviour - clearly visible between 1000 and 
1700 minutes during the 5-day session. For the epoch shown in 
subplot (d), which lasted only for 800 minutes, we were just "un- 
lucky" to catch this kind of variation and not the regular scintil- 
lation pattern. The same applies for the last session, with average 
flux value also exceeding the general average calculated from all 
sessions (see Table[T]). 

Because these last two sessions were so untypical for this 
pulsar, we decided to calculate the general average flux value, 
and also the average modulation index, with and without them. 
Both results are shown in Table[T] This table presents the results 
of the mean flux measurements, averaged over each observing 
session {F}, along with their standard deviations ((T</r>), which 
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Fig. 1. Top left plot: flux time series for the longest observation time (5 days). Top right: flux time series with the largest flux 
variations. Bottom left: one of the unusual flux time series (gaps and some uncharacteristic flux behaviour). Bottom right: strange 
flux variations characteristic for the last two of our sessions. The bar in the top-left corner of each subplot represents Icr flux error 
for each session. 



also include the contribution from the calibration source uncer- 
tainty. We also present the flux modulation indices calculated for 
individual sessions (see EquationlT), these will be discussed in 
detail in the following sections of this paper. 

2.2. Individual session flux density variations. 

The modulation index m is defined as the ratio of the RMS devi- 
ation to the mean value of the observed flux densities 



^|{(F-{F}f) 



(F) 



(1) 



Table [U presents the values of the measured modulation in- 
dices for each individual session. The average values of m in the 
bottom section of the table, just as the flux density averages (see 
previous section), are calculated with and without last two ses- 
sions. For these sessions the quasi-stable flux behaviour means 
that the relative strength of the modulation drops and also causes 
the increase of (F). Both these effects contribute to the decrease 
of the measured modulation index for the two ultimate sessions. 

For the remaining sessions the modulation index varies from 
0.71 (November 2002) to 1.32 (August 2002), with an average 
value of m - 0.96 + 0.04, which means that the observed flux 



density of the source is undergoing significant variations. One 
has to ask a question: what contributes to the observed flux mod- 
ulation? 

We expected that at the frequency of 4.8 GHz 
PSR B0329-I-54 will be still in strong scintillation regime, 
which means that its signal will be undergoing diffractive as 
well as refractive scintillations. This came from the observa- 
tional experience with this pulsar and strong flux variability 
observed at 5 GHz, when we observed it for the purpose of 
estimating its spectrum (Maron et al. 2000). The theory (see 
Lorimer & Kramer. .2005. for summary) predicts that at 4.8 GHz 
a pulsar with DM = 26.8 pc cm"^ should be close to its 
transition frequency, but on the strong scintillation side as well. 

For the DISS the intrinsic modulation index should be close 
to unity. In usual, low-frequency observations, where the decor- 
relation bandwidth Avdiss is significantly less than the observ- 
ing bandwidth, the number of observed scintles within the band- 
width, at any given time, will be large. That would lead to the 
decrease of the modulation of the average pulsar flux (integrated 
over the whole bandwidth), as at any given time, one would see 
multiple scintles, possibly in different stages of development. 

Our observing bandwidth of 192 MHz is at least four times 
more narrow than the expected Avpiss (see Table 3). It means 
that at any given time we are usually able to see only a fraction 
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Table 1. Average flux density measurements for each epoch and 
modulation index (as calculated from Eq[T]i- 



Date 



t obs 

(niin) 



(mJy) 



0"<f> 

(mJy) 



2002-Jul-03-08" 


6300.40 


19.1 


2.0 


1.06 


2002-Aug-22-23'' 


1454.30 


17.4 


2.1 


1.32 


2002-Sep-18-19'' 


1816.80 


27.6 


3.1 


1.13 


2002-Nov- 19-20 


1731.50 


14.0 


1.5 


0.71 


2003-Jan-13-14'' 


1801.30 


18.9 


2.3 


1.28 


2003-Feb-06-07 


1277.00 


20.6 


2.3 


0.83 


2003-Feb-19-21 


2650.30 


11.9 


1.3 


0.77 


2003-Mar-16-18 


1707.50 


6.6 


0.7 


0.88 


2003-May-01-02'' 


1852.30 


11.5 


1.2 


0.91 


2003-Jun-12-13 


1766.00 


13.7 


1.5 


0.75 


2003-Sep-lO-ll 


1921.50 


16.5 


1.8 


1.02 


2003-Oct-21-22 


1684.83 


27.4 


3.0 


0.89 


2004-Jan-03-05 


2248.34 


22.8 


2.4 


0.87 


2004-Mar-15-17'' 


2279.67 


25.8 


2.8 


0.91 


2004-May-ll 


1105.00 


15.7 


1.9 


1.11 


2004-NOV-16-17 


861.00 


6.3 


0.8 


0.93 


2004-Nov-27.11 


519.17 


5.2 


0.7 


1.00 


2005-Feb-16-17* 


781.66 


7.2 


0.8 


0.90 


2005-Mar-26*'' 


801.17 


49.5 


5.4 


0.59 


2005-JU1-24-25'- 


575.33 


30.6 


3.3 


0.53 



Average values: 
all sessions 
w/o last two 



17.4: 
16.0: 



2.4 
1.7 



0.92 ± 0.05 
0.96 ± 0.04 



Notes: " both scintillation timescales observed 

* structure function does not saturate for this 
'^ untypical flux variations - see FigHJl 
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Fig. 2. Flux density for each epoch of flux density monitoring at 
4.8 GHz. The error of the mean flux density was calculated from 
(tF/F X mRiss/(fofc.!/fRiss)'''^ (KS92). The top horizontal Hne 
represents the arithmetic average flux, the lower line corresponds 
to the weighted average flux. 



To take that into account, we calculated the errors of average 
flux measurements following Kaspi & Stinebring ( 119921 KS92) 



~F 



niRiss 



yTobJtRiss 



(2) 



of a single scintle. Following Cordes & Lazio (199^, we esti- 
mated our number of scintles in both frequency and time, and 
calculations yielded Nf to be very close to unity (as expected) 
and Nt - 4.6. The latter leads to the expected DISS modulation 
index of 0.466, much lower than the observed values. However, 
after including the contribution from RISS (wriss = 0.56; the 
value found via structure function analysis, see next sections of 
the paper), and using the total intensity variance formula (Rickett 
I1990I I. we obtained the final expected value of mtot - 0.889. This 
value is still somewhat lower than the observed modulation in- 
dices, but not by a huge margin. 



2.3. Average flux density variations 

Because every individual session was at least several hours long, 
we believe that diffractive scintillations (which happen at the 
timescale of ca. 40 minutes, see subsection 12.4b should not af- 
fect our average flux measurements. On the other hand, refrac- 
tive timescales, which are significantly longer, may have affected 
our results. 

Figure |2] shows the results of the average flux density mea- 
surement versus the observing epoch. Clearly (conf. Table 1) 
there is significant variation in the average values. Using these 
data, we calculated the modulation index of the average flux 
density measurements, which yielded the value of m = 0.57. 
Assuming that the pulsar itself is not varying in intensity, this 
modulation could be caused by only the refractive scintillations, 
which happen at significantly long timescales, close to the length 
of a single observing session. 



where rob,s is the length of a given observing session, and fRiss 
is the RISS timescale. To calculate the error values shown on 
Figure |2] we used the value of fRiss obtained via the structure 
function analysis (305 minutes, see subsection l2.4l i. 

Using those error estimates, we were able to calculate the 
weighted average flux density (w, = cr7^) for our entire observ- 
ing session, which is (Ftot) - 1 1 -69 mJy. This value is shown 
in Figure |2] as a dashed-dotted horizontal line, along with the 
simple arithmetic average from Table [1] (dashed line). 

As we mentioned above, the average flux values and their 
respective uncertainties for the last two sessions in our project 
differ significantly from the remaining observations. However, 
a proper calculation of the eiTor estimates improves the picture. 
For relatively short sessions, which we had towards the end of 
the project, RISS can play a huge role, but the duration of the 
session is included in the uncertainty estimates, which makes 
them more reliable than the formal errors cited in Table 1 . The 
weighted average is therefore much lower than the arithmetic 
average. 

At this point one has to note that the estimation of the flux 
density of pulsars at high observational frequencies will always 
be affected by strong scintillations (except for some very low 
DM cases) and one has to take it into account when trying to 
measure the flux of these, for example for the purpose of acquir- 
ing pulsar spectra. 

2.4. Structure function analysis: timescales 

To estimate the scintillation timescales, we used the structure 
function analysis of the flux density variations. For uniformly 
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sampled data the first order of the normalized structure function 
of the flux density is (Simonetti et al. 1 19851 1 



D(j) 



1 



{FYmj) 



2 w{i)w{i + j)[F(i) - F(i + j)f 



(3) 



where / is the data index, /' is the offset that corresponds to a 
time lag in the data for which the value of the structure function 
is calculated. w{i) is the weighting factor (which is equal to 1 
where the flux density exist for jth points and otherwise), and 
N(j) = 2 w(i)w(i + j) is the number of pairs of values found in 
the data that were found to have the offset equal to j (i.e. time 
lag equal to r). 

The values of the structure function need to be corrected to 
remove white noise contribution. This was made by subtracting 
the value of the structure function at the unit lag from all the 
structure function values. If a quasi-periodic signal is present in 
the data the structure function will show a plateau at the value 
of Dsat. The timescale of the observed variability is the lag cor- 
responding to half the saturation value of the structure function 
after white noise correction: fiss = T(Dsat/2) (KS92). 

The error of the timescale can be estimated by finding the 
time lags that correspond to the values of the structure function 
of (I>sat ± <5Dsat) /2. The uncertainty of the saturation value of 
the structure function SD^ai is best estimated as SDsailDsat ~ 
(2fiss/robs)'^^, where Tobs is the total duration of the observa- 
tion, in our case the length of an individual observing session. 

Since our project consisted of several observing sessions re- 
peated over the course of three years, we decided to calculate 
the general average structure function as well. This can be made 
either directly by applying structure function algorithm to the 
entire set of data, or by adding the values obtained for individual 
sessions, after de-normalizing them. 

The values of the scintillation timescale, modulation index 
and their uncertainties for the general average structure function 
can then be calculated in the same manner as for the individual 
sessions. 

Table |2] shows the results of the structure function analysis 
applied to our data. This table has less individual session en- 
tries than Table 1 because it shows only those sessions for which 
the structure function saturated. The saturation of the SF was 
necessary to calculate the cited values of fiss and m and their 
uncertainties (see the notes for Table 1). 

Figure |3] shows three of the structure functions obtained dur- 
ing the analysis. The top plot shows the SF for the first and 
longest, 5-day session that was conducted between 2002 July 
3-8. This structure function clearly shows two plateaus, which 
allowed us to obtain the values of two distinctive timescales 
present in the data. As we mentioned above, at the observing fre- 
quency of 4.8 GHz PSR B0329+54 is beheved to be close to the 
transition frequency, i.e. switching from strong to weak scintilla- 
tion regimes, difFractive and refractive timescales should be rela- 
tively close (yielding a low value of the strength of the scattering 
parameter u). Hence we are convinced that those two plateaus 
correspond to two scintillation timescales, a shorter diffractive 
timescale (lower plateau), and a longer refractive timescale. 

This is even more convincing for the middle plot of 
Fig. [3] which shows the general average structure function. Two 
plateaus were also present for another session, 2002 September 
18-19 (ca. 30-hour session). 

The third and bottom subplot of Fig. |3]shows a typical struc- 
ture function, obtained during a single, typical 37-hour session. 



Table 2. Parameters calculated from structure function for flux 
density series: scintillation timescales and flux modulation in- 
dices. The top part of the table presents the values for DISS 
timescales and modulation indices for individual sessions, as 
well as mean values and results from general average SF analy- 
sis. The bottom part presents RISS parameters for two individual 
sessions and the general average SF. 



Date 


Lag 


?ISS 


cr(fiss) 


m 


0-,,, 




(min) 


(min) 


(min) 






2002-JuI-03-08 


6300 


58.3 


+5 

-A 
13 


0.95 


0.15 


2002-Aug-22-23 


1454 


45.0 


1.23 


0.34 


2002-Sep-18-19 


1816 


29.0 


+4 


0.89 


0.2 


2002-Nov- 19-20 


1731 


33.4 


+5 
10 


0.70 


0.2 


2003-Feb-06-07 


1227 


45.8 


0.81 


0.29 


2003-Feb- 19-21 


2650 


33.7 


+6 
''O 


0.77 


0.16 


2003-Mar-16-18 


1707 


32.4 


0.82 


0.21 


2003-Jun-12-13 


861 


29.5 


0.80 


0.27 


2003-Sep-lO-ll 


1921 


58.1 


1.00 


0.29 


2003-Oct-21-22 


1684 


81.6 


0.91 


0.35 


2004-Jan-03-05 


1489 


34.5 


+3:5 

10 


0.86 


0.23 


2004-May-ll 


734 


37.0 


0.90 


0.36 


2004-NOV-16-17 


861 


37.4 


+7 
-4 


0.93 


0.34 


2004-NOV-27 


519 


25.4 


1.02 


0.38 


2005-Jul-24-25 


575 


19.0 
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where one can see only the regular saturation of the structure 
function caused by the difFractive scintillations. 

The horizontal lines in all those plots show the SF saturation 
levels, vertical lines - the time lags that correspond to the satu- 
ration (the values of the timescales are given in minutes), while 
the short horizontal dashes crossing them represent the errors in 
the timescale estimates. 

All detailed results are presented in Table |2l The values of 
fRiss for individual sessions differ from the average SF values. 
This is understandable because the RISS saturation level is de- 
termined with lower precision, because it is limited by the max- 
imum time lag of the SF calculated for the given session, and 
because especially in the September 2002 we saw only 4-5 RISS 
cycles. This translates into a large uncertainty of fRiss. While for 
the September 2002 session the values agree with those obtained 
from the general average SF within error estimates, for the first 
and longest July 2002 session they do not. This may be because 
this session may also suffer from the untypical flux variations 
that were observed (see Figure 1 and Section lTTT i. 

Because the general average structure function of our data 
is the best representative for our results, we'll be using those 
values for further calculations: foiss - 42.7 minutes, fRiss = 
305 minutes. Given these, we can calculate the strength of the 
scattering parameter u - V^Riss/foiss - 2.(fl . 

Also, since in our observations we were unable to obtain the 
value of the decoiTelation bandwidth Bdiss (this value greatly 
exceeds the total observing bandwidth of PSPM 2, which is 
192 MHz) we used the results of our timescale measurements 
to estimate this value. Using the Stinebring & Condon (II990I I 
formula: 



Bdiss - 1-27/obs 



^DISS 
fRISS 



(4) 
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Fig. 3. Three structure functions of our data. From top to bot- 
tom: the structure function from the 5-day session, general av- 
erage structure function, and the typical structure function of an 
individual session (see text for details). 



we obtained a value of Z?diss - 853 MHz (over 4 times higher 
than the observing bandwidth). 



2.5. Modulation indices from the structure function analysis 

As we mentioned in the previous section, the structure func- 
tion analysis allows us to calculate not only the scintillation 
timescales, but also their respective scintillation indices. For 
a normalized and white-noise corrected structure function the 
value of the scintillation index m = (Dsat/2)''^, and the frac- 
tional error in m is a sum of errors from the fractional errors in 
Dsat and m (see KS92). 

Table |2] shows the results for the individual sessions as well 
as for the general average SF for the modulation indices of both 
difFractive and refractive scintillations. 

For the difFractive scintillations the average value of moiss, 
obtained from the simple arithmetic average of the individual 
values agrees with the value obtained from general average 
structure function, which is equal to wjdiss - 0.89. This value 
differs slightly from the modulation index value obtained di- 
rectly from flux density measurements (0.96, see Table 1), but 
is very close to the expected value for the observed modulation 
index (0.889, see section |2!2] |. 

The values of the RISS modulation indices are also presented 
in Table|2] For individual sessions that allowed for an estimation 
of wjdiss, they diff'er slightly from the average values, but still 
agree within the error estimates. This is not a surprise, because 
it is diflicult to determine the saturation level properly, especially 
for the shorter sessions (for the same reasons the timescales 
differ, see previous section). On the other hand, the value of 
wiRiss = 0.56 obtained for the general average structure func- 
tion is very close to the value of the RISS modulation index ob- 
tained from the average flux density variation analysis (0.57; see 
section |23] and Fig.|2]i. 

Again, as in the case of the ISS timescales, we will be using 
the values obtained for the general average structure function as 
the final result of our analysis: moiss = 0.89 and wriss - 0.56. 



3. Discussion 

The results obtained from the analysis of our observations sug- 
gest that we definitely see both refractive and difftactive scintil- 
lations affecting the measured pulsar flux density. We were able 
to calculate both the timescales and modulation indices for these, 
and also to estimate the expected decorrelation bandwidth. 

Our observations were performed at the frequency of 
4.8 GHz, for which there is only one observational report to be 
found in the literature by Malofeev at al. ( 1996). This was quite 
a short project, because PSR B0329-I-54 was only observed for 
130 minutes at this frequency, which given the values of the scin- 
tillation timescales we found makes the results unreliable. On 
the other hand, adding our results to the scintillation parameters 
estimates obtained by other authors at lower frequencies allows 
us to study the frequency dependence of these parameters over a 
much wider range of frequencies. 

3. 1 . Flux density measurements 

As we mentioned above, PSR B0329-I-54 at the observing fre- 
quency of 4.8 GHz shows large flux density variations owing to 
ISS. For some of the individual 3 -minute integrations the flux 
measurements yielded values as high as 200 mJy, over 10 times 
higher than the average value. Those variations happen at a va- 
riety of timescales, ranging from minutes to hours, and up to 
days (c.f. Fig la). This clearly shows the dangers of estimating 
pulsar flux densities at high frequencies (especially for low-DM 
pulsars) from observations that are based on a limited number 
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of separate integrations, which is often the case when one tries, 
for example, to ascertain the pulsar spectrum. The diffractive 
timescale is longer at high frequency, and may be comparable 
to the actual integration times used - usually of the order of a 
few tens of minutes when one tries to measure the flux density 
of some weaker sources - hence, it may strongly affect the out- 
come. On the other hand, with refractive timescales of the order 
of a few hours, repeating the observations after a day or a few 
days, with the hope that the ISS will be averaged-out, may not 
solve the problem. One has to approach these cases carefully; 
luckily, for high-DM pulsars this should be less of a problem, 
because the transition frequency will be much higher and the 
modulation caused by RISS should be very low, and one has to 
worry only about diffractive scintillations when performing flux 
measurements. 

In case of our observations of PSR B0329+54 we were 
able to ascertain the impact of interstellar scintillations on 
our measurement pretty well. Both our values: 18.4 mJy 
from the simple arithmetic average and 11.7 mJy from the 
RISS weighted average agree well with the simple power- 
law spectra for this source, with spectral index of -2.2 (see 
Maron et al. 2000, and the authors on-line spectra database at 
http://astro.ia.uz.zgora.pl/olaf/paperl/index.htmI). 

3.2. Scintillation parameters 

Pulsar B0329-H54 is one of the strongest radio pulsars known 
so far, and scintillations of its radiation were analysed at various 
frequencies by many different authors, especially over the course 
of last 20 years. Tableland Figure |4] summarize the results we 
were able to find. 

Table [3] shows the list of all scintillation parameters for 
PSR B0329-I-54 at diff'erent frequencies ranging from 74 MHz 
to 4.8 GHz. At 4.8 GHz, at which we conducted the ob- 
servations as well, we only found one measurement of foiss 
in Malofeev et al. (1 19961 1. which is 21.4 minutes, almost ex- 
actly two times lower than our value from the general aver- 
age structure function. This discrepancy may arrise because 
Malofeev et al. (1996) observations were very short: 130 min- 
utes only, i.e. three diffractive cycles, according to our re- 
sults. Additionally, one has to note that even during our ob- 
servations we found that for individual sessions the measured 
value of fpiss can go as low as 19.0 minutes (see Table |2]). 
Malofeev et al. (1996) also estimated a lower limit for the 
diffractive modulation index (moiss > 0.5). 

Figure |4] shows the frequency dependence for both diffrac- 
tive and refractive scintillation timescales (top right and top left, 
respectively). Data taken from the literature are represented by 
full circles (for reference see Table [3]), our measurements are 
the open circles. Evidently, our measurements for both foiss and 
fRiss roughly agree with the previous, lower frequency estimates 
in a sense that they support the theory that the frequency de- 
pendence of the scintillation timescales is really governed by a 
simple power-law. Adding our points to these plots extends the 
frequency range significantly, which in turn allows for better fit 
of the power-law to the data. This is especially important for 
the RISS timescale, because it is very long at lower frequencies 
and the measurements often suffer because they are either poorly 
sampled or have only short time spans; i.e. do not cover a suffi- 
ciently large number of observed refractive cycles. These result 
in a large spread of fRiss measurements at lower frequencies. 

On the other hand, our measurements were obtained from 
the continuous observations that usually covered at least 4-5 re- 
fractive cycles (in a normal one-day session) and up to over 20 



cycles (for the longest 5-day session), which makes them very 
reliable. 

Using the data summarized in Table [3] we performed power- 
law fits, which yielded foiss = 9.3 fj,f^ (minutes) and fRiss = 
3.9/^,-1™ (days). 

Similar fits were performed for the RISS modulation in- 
dex (see bottom left plot in Fig. |4]) and the decorrelation band- 
width BdisS' although the latter is not a direct measurement 
but only an estimate (see Section |24] ). Nevertheless, both mea- 
surements roughly agree with expectations based on lower fre- 
quency data, and adding our points to the pool allowed us to 
obtain the frequency dependencies as niRiss - 0.38 fj^'^^ and 
Bdiss = 2.6 f2f. 

We have to note that for the purpose of two of these fits, 
namely Jriss and wriss, we decided to omit the values obtained 
by Wang et al. ( 120081 ). which are represented by triangles in their 
respective plots. These points stray by much from the rest of 
the data, and contradict the earlier estimates at the same fre- 
quency (1540 MHz) given by Wang et al. (I2005I) . Adding our 
values to the data pool clearly shows that their earlier results (ob- 
tained from a different set of observations) seem to fit better into 
the general picture. We believe that the 2008 results are flawed; 
problem may lie in the method that was used for the purpose of 
that paper. Based on other data, the actual refractive timescale 
at 1540 MHz should be of the order of 2 days (similar value as 
that given by their earlier paper, which is 2.5 days). Their ob- 
servations were conducted almost continuously for 20 days in 
March 2004, with the flux measurements based on 90-minute 
integrations. At the frequency of 1540 MHz this means that they 
observed only a few diffractive cycles during their single flux- 
integration, which may have affected the measurements. The au- 
thors calculated wriss based on that data. However, their struc- 
ture function for the flux measurements did not saturate. Without 
the saturation level it is impossible to estimate the timescale 
from the structure function (see section 12.41 ). and they decided 
to obtain the saturation level as Dsat - in^/2 (see section 1231 ). 
Using the value obtained this way, and applying it to the struc- 
ture function (see their Fig. 5), they inferred the value of the 
RISS timescale of 8 hours, 6 times lower than expected and al- 
most 8 times lower than in their previous paper. 

Probably the flaw lies in the underestimation of the RISS 
modulation index (which clearly shows in our Fig. |4] as well) 
which resulted in a lowering of the used saturation level. An in- 
spection of the structure function they obtained (see their Fig 5, 
Wang et al. 12008!) leads to the conclusion that even a slight 
change of toriss, and in turn the saturation level Dsat, could lead 
to timescale even a half an order of magnitude higher than the 8 
hours they quote, and we think that is exactly what happened. In 
Wang et al. (120081) defense we may note that we had one session 
in our data that is coincident with their March 2004 observations, 
as we performed a 38-hour observation between March 15 and 
17 of 2004, and our structure function did not saturate as well, 
which may suggest that the fluctuations of the flux density were 
untypical during that period. 

Knowledge of frequency dependence of both foiss and Jriss 
gives us an opportunity to calculate the transition frequency, 
i.e. the observing frequency at which both timescales would be 
equal, which means a switch from strong to weak scintillation 
regimes. Using the empirical formulae provided above (see also 
FiglH and putting both timescales equal, this yields the transition 
frequency Vc of 10. 1 GHz. Below this frequency PSR B0329-H54 
will be in the strong scintillation regime, above that frequency. 
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Fig. 4. Frequency dependence of diffractive timescales foiss (top left), refractive timescales tjass (top right), refractive modulation 
index toriss (bottom left), and decorrelation bandwidth BW wriss (bottom right). Data taken from Bhat et al. dl 9991 1 at 327 MHz, 
Gupta et al. ( 1994) at 408 MHz, Stinebring et al. (fT996l l at 610 MHz, Wang et al. (.2005J at 1.5 GHz, and this paper at 4.8 GHz. The 
sohd hne presents the fit to observed data. 



the characteristic of the scintillation should change to weak scin- 
tillations. 



3.3. Estimation of electron density turbulence spectrum 

Observations of the interstellar scintillations and especially the 
frequency dependence of the measured parameters can be used 
to estimate the properties of the interstellar turbulence. By anal- 
ogy to the neutral gas theory, the density fluctuations in the ISM 
can be described by a power-law spectrum as P^n - Cl q^^, 
where C^ is the mean turbulence of electron density along the 
line of sight, q - Inj s is the wavenumber associated with the 
spatial scale of turbulence s, and the spectral index fi, which 
is in range of 3 < /? < 5. The well known Kolmogorov the- 
ory developed for neutral gas turbulence yields the spectral in- 
dex of yS = 11/3, but it was shown by several authors (see 
for example Gupta et al. '1993) that there are discrepancies be- 
tween Kolmogorov theoretical predictions and actual measure- 
ments of scintillation parameters for some pulsars. Few mod- 
els of the interstellar turbulence with different spectral slopes 
P were developed, yielding different sets of predictions of the 
frequency dependence of ISS parameters (see Romani et al. 
[1986 ; Bhat et al. 1999). Table |4] summarizes these predictions 
for the three most commonly used spectral slopes of y6 = 1 1/3 
(Kolmogorov spectrum), 4 and 4.3 for the four scintillation pa- 



rameters we were able to measure, or derive from our observa- 
tions of PSR B0329+54. Our estimations of the actual spectral 
slopes of these parameters, based on all available data (see sec- 
tion l3.2l and FigurelH are also in the table. 

Apparently none of the currently models for the electron tur- 
bulence agrees precisely with our estimations, and the model 
predictions with (3 - A are the closest to the slopes we ob- 
tained from our fits. It has to be noted that in case of three 
of the four parameters (with the wjriss being the only excep- 
tion), these discrepancies cannot be attributed to poor quality 
of frequency slopes fits for these parameters, especially that our 
measurements, made at the frequency of 4.8 GHz, significantly 
widened the frequency range, putting strong constraints on the 
fits. 



3.4. Derived scintillation parameters 

As we mentioned above, in a strong scintillation regime the flux 
density variations happen at two distinctive timescales - diffrac- 
tive and refractive. The apparent variability of the pulsar sig- 
nal can be understood as the result of the observer crossing the 
diffraction pattern, which is introduced by the ISM to the pulsar 
signal wavefront. In a thin screen model, those two timescales 
correspond to two spatial scales of the diffraction pattern: the 
diffractive scale s^i, and refractive scale Sr, which in turn can 
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Table 3. Estimated and observed parameters. Numbers in parentheses are error estimates at the level of the last digit(s) quoted. 



foiss 
(min) 


?RISS 

(min) 


"JDISS 


'MRISS 


u 


Bdiss 
(MHz) 


VlSS 

(km/s) 


References 
and notes 


Our 4.8 GHz Results 42.7 

derived parameters 


305 


0.87 


0.56 


2.67 


853 


92.2 




Results from literature: 



4.8 GHz 
1540 MHz 
1540 MHz 
750 MHz 
610 MHz 
416 MHz 
408 MHz 
327 MHz 
327 MHz 
310 MHz 
74 MHz 



21(4) 

16.9 
6.64 
5.9 

3.27 
3.23 
4.57 

2.3 



(?)480 
3600 

4.0 days 
15 days 

14.7 days 



31(14) days 

28 days 
> 100 days 



>0.5 



0.47 



(?)0.31 



0.37(3) 



0.27(7) 
0.15(6) 



15 



63 



14.0 
1.096 
0.349 
0.081 
0.047 
0.165(8) 

0.0022 



1 (short observations) 

2 (doubtful) 
97 3 

4 

5,6 

4 

7 

8 

9 

4 

7 



References: (1) Malofeev et al. n996^, (2) Wang et al. i BOOSl l. ( 3) 154 MHz Wang et al. f2005^, (4) Stibenring & Condon l fT990t . 
(5) Stinebring et al. Cl996j . (6) Stinebring et al. ( l2000t . (7) Gupta et al. ( fT994> . (8) Bhat Cl999j , (9) Esamdin et al. (.2004J 



Table 4. Theoretical and observed spectral indices of scintilla- 
tion parameters (Romani et al. .l986. Bhat et al. 1 19991 1 
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■al index predicted 
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by theory with yS 
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spectral index 




11/3 




4 




4.3 


(our work) 


'diss 


1.2 




1.0 




1.4 


1.01 


fRISS 


-2.2 




-2.0 




-2.4 


-1.76 


"JRISS 


0.57 




0.38 




0.55 


0.32 


Sdiss 


4.4 




4.0 




4.7 


3.68 



be bound to the concepts of scattering angle and scattering disk 
that is commonly used in the scintillation theory (see Rickkett, 
fT990l l. 

Following Gupta et al. jl994l l. we can find the diffractive 
scale by the means of estimating the scintillation velocity Viss, 
which is bound to the diffractive scale as s^ - V/ssfoiss- 
The scintillation velocity can be derived from the diffractive 
timescale fniss and the decorrelation bandwidth Bdiss (Cordes, 
[1986', Gupta et al. 1 19941 ). Using our values in the formula yields 
Viss = 92.9+^0^4 km/s. 

The observed scintillation velocity should be modulated by 
the Earth's orbital motion, but our data proved to be inconclusive 
in that regard. We only had the diffractive timescale measured for 
each epoch, and were forced to use the average value of decorre- 
lation bandwith (and/or refractive timescale) in the Viss formula, 
which led to a wide spread in the Vjss versus epoch data, making 
our fits unreliable. Measurement of the refractive timescale for 
every epoch would probably solve the problem, as would the di- 
rect measurement of the decoiTelation bandwidth, but these were 
impossible owing to the character of the variability of the pulsar 
flux and the limited observing bandwidth of the backend used. 

The average value of Viss - 92.9 km/s combined with the 
measured foiss = 42.7 min. yielded the diffractive scale Sd - 
2.38 X 10^ meters. 



The diffractive 
(Wang et al. l2005T l 



Od 



n d Buiss 



1/2 



angle 6^ can be estimated from 



(5) 



where c is the speed of light. This yields 0^ = 0.01207 milliarc- 
seconds, which translates into the refractive scale Sr - 1 .91 x 10'^ 



meters. Combining this with the previously obtained value of 
diffractive scale sj, we can estimate the Fresnel scale for the ob- 
served ISS as rf - (sds,y^^ - 6.7 x 10^ meters. 

One has to remember that for observations we were not able 
to measure the decorrelation bandwidth Bdiss directly, and both 
of the formulae for Vjss and 0^ are based on this value. 

The value of the strength of the scattering parameter inferred 
from the spatial scales s,- and Sd i& u = 2.84, which is slightly 
different from the value derived from the scintillation timescales 
(u - 2.67, see Section lZ^l l. The discrepancy is small, however, 
and given the fact that the application of the ISS theory to the ob- 
servations always involves lots of approximations and assump- 
tions, our values of Sr, Sd and rf we got are as reliable as possi- 
ble. 

4. Summary and Conclusions 

We presented the results of our three-year observing project of 
PSR B0329+54 at the frequency of 4.8 GHz, using the 32-meter 
Toruii Centre for Astronomy radiotelescope and the Penn State 
Pulsar Machine II as a backend. The project was very unique, be- 
cause it consisted of twenty separate observing sessions, each of 
them involving several hours (up to 5 days; more than 24 hours 
on average) of continuous observations of the flux density of 
the pulsar, with 3-minute time resolution. The total observing 
time ammounted to 35 thousand minutes, or 24 days, making it 
probably one of the longest observing programmes that involved 
long-term continuous observations of flux density and scintil- 
lation parameters that was performed using a single observing 
setup and uniformly analysed from raw data to the final results. 

4.1. Scintillation parameters at 4.8 GHz 

The analysis of the data allowed us to conclude that flux density 
variations are governed by two distinctive timescales, which we 
identify as the diffractive and refractive scintillation timescales. 
The fact that the observing frequency of 4.8 GHz is relatively 
close to the transition frequency for this pulsar makes those 
timescales differ only by a factor of ~ 7. The character of our 
observations allowed us to detect them both at the same time, 
by means of the structure function analysis, using the same set 
of data - twice for individual sessions and additionally by con- 
structing general average structure function. To our knowledge 
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this is the first case of pulsar ISS observations where both scin- 
tillation timescales were measured from a single data set using 
exactly the same method. 

The measured values of the diffractive timescale foiss - 
minutes and refractive timescale fRiss - 305 ^^^^ minutes 
allowed us to estimate the decoiTelation bandwidth (unobtain- 
able via direct measurements, because it exceeds the observing 
bandwidth of the backend used) as Bdiss - 853 MHz. Also, 
we were able to estimate the strength of scattering parameter, 
which can be expressed as square-root of the timescale ratio, as 
u = 2.67. 

We were also able to estimate the amount of modulation that 
both types of scintillation contribute to the flux variability, and 
did that by the means of the modulation indices. Especially im- 
portant is the estimation of the modulation due to refractive scin- 
tillations, for which the modulation index is relatively strong: 
'Wriss = 0.56. This agrees with the fact that during some of our 
observing sessions we noticed very drastic flux density varia- 
tions, with the flux from single 3-minute integration reaching 
the values of 200 mJy, exceeding the average value by a factor 
of ~ 20. 

This shows how important it is to take ISS into account 
when measuring pulsar flux densities at high frequencies, espe- 
cially for low-DM pulsars. For these objects the frequency of the 
switch from strong to weak scintillation regimes (i.e. the transi- 
tion frequency) should be relatively low, and all theories pre- 
dict the increase of refractive modulation index towards that fre- 
quency (for summary see Lorimer & Kramer 2005 ; also Romani 
et al. I1986I I. In these cases, when the ISS timescales are of the 
order of tens to a few hundred minutes, typical flux density ob- 
servations (which usually involve continuous integration of the 
length of 10 minutes to an hour) of a single measurement may be 
heavily affected by both types of scintillations. This in turn may 
cause problems when trying to construct the pulsar spectra, lead- 
ing to artificial breaks or turn-ups observed in the spectrum - as 
was pointed out in a few cases by Maron et al. (2000). The only 
way around the problem is to repeat the measurement several 
times to average-out the influence of ISS. Luckily, this should 
be less of a problem for pulsars with higher dispersion mea- 
sures because the transition frequency for such objects should 
be well outside the usual observing frequencies, meaning that at 
least the refractive scintillations should not affect the measure- 
ment by much. Still, the diffractive scintillations would present 
a problem if their timescale would be similar to the integration 
lengths used, so one should approach every case carefully and 
individually. 

Our measurements of the scintillation timescales allowed 
us also to estimate some of the parameters commonly used to 
describe the theory of interstellar scintillation. We were able 
to estimate the Fresnell scale rp as 6.7 x 10^ meter, which at 
the frequency of 4.8 GHz corresponds to a refractive scale of 
Sr - 1.9x10^ m, and the diffractive timescale of sj - 2.4x10'^ m. 

We were also able to estimate the scintillation velocity of 
Viss - 92.9 km/s, and we hope that our value will help to solve 
the problem with the discrepancies that can be found in the lit- 
erature. It agrees well with Wang et al. ( I2005I I as they obtained 
97 km/s, on the other hand KS92 report Vjss - 68 km/s, whereas 
Gupta (119951 1 170 km/s or 126 km/s, depending on the model. 
One has to remember that those discrepancies may arise because 
while the the formula is uniform (with the possible exception of 
the scaling factor A,,), the ways used to obtain the parameters 
were usually different. 



4.2. General picture 

Adding our results to the estimations of the scintillation param- 
eters at lower frequencies that can be found in the literature al- 
lowed us to study the dependence of these parameters on the 
observing frequency. That our observations were conducted at 
the frequency of 4.8 GHz proved to be very useful because it 
significantly widened the frequency range for which the values 
of ISS parameters are known. Previous attempts at this analysis 
were limited to the frequency of ~ 1.5 GHz, adding our mea- 
surement expanded that range three-fold, which means an exten- 
sion of ca. 50% when talking about the log-scale in frequency. 
This proved to be especially helpfull for the RISS parameters. At 
low frequencies the RISS timescale is very long and the modu- 
lation very weak, which brings a relatively large scatter to the 
low-frequency parameter estimates, making any attempts to fit 
a power-law type dependency questionable. The addition of our 
results, which we believe to be very reliable (because they come 
from a long-term monitoring programme), improves the general 
picture greatly (see FigHJl, making any fits to the data much more 
constrained. 

We performed the fits of the power-laws governing the fre- 
quency dependence of the scintillation parameters, and they 
yielded us the respective spectral indices, which can be ex- 
pressed as: fDiss ~ / ■•°', fRiss ~ Z"'™, '«riss ~ /"'^^ and 
Bdiss ~ / ^^^- These values significantly differ from those pre- 
dicted by using Kolmogorov neutral gas theory to the ISM spa- 
tial electron density spectrum (fi - 11/3, see Table |4|i. Indeed 
our measurements do not agree with any of the models de- 
scribing the electron density spectrum that are commonly used, 
with the model using fi - A being the closest to the obtained 
spectral indices. This would indicate that the real spatial elec- 
tron density spectrum is steeper than the Kolmogorov spectrum. 
Unfortunately, the character of our data did not allow for the es- 
timation of the y6 index by the means used by a few other authors 
(see for example Wang et al. 2005) because they usually involve 
measurements performed in the dynamic spectra. 

Using the above mentioned formulae we were able to esti- 
mate the transition frequency for PSR B0329-H54, i.e. the fre- 
quency at which the transition from strong to weak scintillation 
mode should appear, which is v^ - 10.1 GHz. This clearly con- 
tradicts the earlier statement of Malofeev et al. ( 119961 ), which es- 
timated the transition frequency to be around 3 GHz. This value 
is also higher than the transition frequency one can infer from 
Lorimer & Kramer (120051) - namely from their Figure 4.3, which 
suggests that for a pulsar with DM - 27 pc/cm^ Vc should be 
close to 4 GHz. 

It would be extremely interesting, but on the other hand very 
hard, to perform an observing project similar to ours at a fre- 
quency close to this transition frequency. First - the pulsar would 
be much weaker, thus requiring longer integration times. On the 
strong-scintillation-regime side of the transition frequency, the 
two ISS timescales would be extremely similar (100 minutes), 
which means that the distinction between them would be very 
hard to make, and would probably require very long sets of con- 
tinuous observations. Because there are only a handful of tele- 
scopes that could perform these observations (i.e. have the re- 
quired receiver equipment) this may be impossible to accom- 
plish. 

Nevertheless, any type of scintillation or flux density varia- 
tion observations, performed at the frequency close to the tran- 
sition frequency, may be very helpful to our understanding of 
the scintillation phenomenon in general, and the switch from 
strong- to weak-scintillation in particular This is also true for 
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other pulsars, especially those with lower dispersion measures, 
because the expected transition frequencies for them would be 
even lower than the 10 GHz for PSR B0329+54. 
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